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Abstract 

We develop an iterative approach to gravitational lensing theory based on 
approximate solutions of the null geodesic equations. The approach can be 
employed in any space-time which is "close" to a space-time in which the null 
geodesic equations can be completely integrated, such as Minkowski space- 
time, Robertson- Walker cosmologies, or Schwarzschild-Kerr geometries. To 
illustrate the method, we construct the iterative gravitational lens equations 
and time of arrival equation for a single Schwarzschild lens. This example 
motivates a discussion of the relationship between the iterative approach, 
the standard thin lens formulation, and an exact formulation of gravitational 
lensing. 

I. INTRODUCTION 

In several recent papers an exact approach to gravitational lensing has been de- 

veloped in which a parametric represention of the past light-cone of an observer is viewed 
as the fundamental gravitational lens equations. In this approach, one would solve, in 
principle, the null geodesic equations of general relativity in an arbitrary space-time and 
construct the past light-cone of an observer. The lens and time of arrival equations then fol- 
low from a particular parametric representation of the past light cone. The exact approach 
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can be implemented in only a few space-times possessing a high degree of symmetry, such 
as Minkowski space, Robertson- Walker cosmologies, or Schwarzschild/Kerr geometries. We 
will refer to space-times in which the null geodesic equations can be integrated exactly as 
integrable geometries. 

The exact approach stands in contrast to the generally accepted thin lens approximation 
in which one does not attempt to find the exact null geodesies. In this approximation, 
one introduces sharp bending at isolated points near the lens along a geodesic trajectory of 
a background space-time - usually either Minkowski space or a cosmological model. The 
thin lens approximation has proved extremely valuable to practicing astrophysicists as a 
relatively accurate method which is easy to implement [^. 

Our goal in this paper is to propose a middle ground in which one systematically finds 
approximate null geodesies in geometries "close" to integrable geometries and uses these 
approximate geodesies as the basis for gravitational lensing calculations. These geometries 
need not correspond to any solution of the Einstein equations; our method should apply to 
any approximate metric which can be considered as a perturbation of an exact integrable 
geometry. A key component of this project is a perturbation method known as the "variation 
of constants,' which can be found in several classical mechanics textbooks This method 
is based on techniques from Hamiltonian mechanics and is described in Section II. Though 
it is often used for perturbations of periodic or quasi-periodic orbits, it is suited to finding 
approximations to arbitrary trajectories, regardless of their type. 

In Section III, the method is specialized to the problem of finding approximate null 
geodesies in space-times close to integrable ones, and turning these geodesies into lens 
equations. As an example of the method, we apply the iterative approach to a single 
Schwarzschild lens in Section IV. In a final section, we discuss the relationship between the 
iterative, thin lens, and exact formulations of lensing theory. 
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II. FINDING APPROXIMATE TRAJECTORIES 



In this section, we describe a method to obtain approximate trajectories for a perturbed 
Hamiltonian. The discussion largely follows a similar discussion in Goldstein's Classical 
Mechanics 0. 

Throughout this paper, we refer to an unperturbed Hamiltonian, Hg, and a perturbed 
Hamiltonian, 

= + AH{x^,pa). (1) 

In addition, we will make no distinction between the time coordinate, t, and the spatial 
coordinates, x*. The base space coordinates, x", are thus x"" = (t,a;*), and the momentum 
coordinates, Pa, include Pq, which is canonically conjugate to t. A parameter. A, will be an 
affine parameterization of the trajectory, and is a linear function of t if the Hamiltonian does 
not depend on time. A dot derivative will always refer to a total derivative with respect to 
A, while partial A derivatives will always be explicitly written. 

The method begins by assuming that one can find a complete solution, Fo(x", Pa, A), to 
the Hamilton- Jacobi equation for the unperturbed Hamiltonian, 

with n constants of integration, P^. When 

Fo{x°', Pa, A) is the generator of a (parameter dependent) canonical transformation, 
{x"',Pa) =^ {X"',Pa). The transformation is defined by 

OF 

X- = ^^^X\x\Pa,\) (4) 
dF 

p, = ^^p,(x^P„A), (5) 
with the inversion of Eq.(^) being 



^a^^a^^a^p^^^y (6) 

When the inversion is substituted into Eq.(^, one has the canonical coordinate transforma- 
tion in the form 

Pa = Pa Pa, A), Pa, A) = 7^„(X^ Pa, A) . (7) 

It is important to note that Eq.(0) represents the solution to the equations of motion of 
the unperturbed Hamiltonian in terms of 2n constants, {X"", Pa), which represent the initial 
location and momentum, and the affine parameter, A. 

Since the property of being a canonical transformation is independent of the particular 
form of the Hamiltonian, we may apply the canonical transformation generated by the solu- 
tion to the unperturbed Hamilton- Jacobi equation to the perturbed Hamiltonian. Under any 
canonical transformation generated by the complete solution, -Fo(x", Pa, A), the Hamiltonian 
becomes 

dF 

H\X\ Pa, \)=H {X\X\ Pa, \),Pa{X\ Pa, A)) + (8) 

Using H = Ho + AH, we have 

i^'(X^ Pa, A) = AH (X'^(X«, Pa, \),Pa{X\ Pa, A)) , (9) 

because Foi^x"", Pa, A) was chosen to satisfy the Hamilton- Jacobi equation of the unperturbed 
Hamiltonian. 

Hamilton's equations of motion for X" and Pa are 

r)H' 

X«(X^P„A)= ^(X»,P„A) 
r)H' 

Pa{X\ Pa, A) = -^(^^ Pa. A). (10) 

These equations are exact. We note that if the perturbation, AH, was zero, the solution to 
Eqs.([TO|) would be that X" and Pa are constants. When AH is not zero, these "constants" 



become functions of initial conditions, {X^, P°), and the parameter, A. For this reason, this 
method is referred to as the "variation of constants." 

In principle, the solution of the equations of motion, Eqs.(|TO|), would be 2n functions of 
the form 

x^ = x%x:,p:,x) 

Pa = Pa{XlPlX), (11) 

where and P° are 2n initial values. Substituting Eqs.(|ll]) into the canonical transfor- 
mation, Eqs.(l^), yields the solution to Hamilton's equations for and pa as functions of 
initial conditions and A. 

Note that we have made no assumption about the value of the unperturbed Hamiltonian. 
As an example, suppose that we wish to solve a Kepler type problem, where Ho represents 
a spherically symmetric Hamiltonian and /S.H represents a perturbation. Since the value of 
Ho has not been determined, this method applies to perturbations of bounded or unbounded 
orbits. 

Since the equations for X" and Pa, Eqs.(|TU]), are exact, they tend to be as difficult to 
solve as Hamilton's equations in the original phase space variables, (x",pa). Typically, the 
variation of constants method is not used to find exact solutions, but rather it leads to an 
approximation procedure. 

The natural approximation procedure associated with the variation of constants method 
is to solve Hamilton's equations, Eqs.([10|), by successive iteration. To iterate, one must 
choose an appropriate initial trajectory, (X^(A), P°(A)), as the "zeroth iterate" and insert 
these functions of A into the right hand side of Eq.(|l^). In this way, the right hand sides 
of Eqs.(p!0|) become functions of the parameter A. The first iterate approximation is then 
obtained by simple integration on A. We want to emphasize that any trajectory, suitably 
close to the exact result, can be used as the zeroth iterate. In practice, general knowledge 
of the physical situation will guide the choice of the zeroth iterate. 

One natural choice is to take the zeroth iterate as given by X"- and Pa as constants, 
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X"- = and Pa = Pa, which is the result obtained when AH = 0. The first iterate then 
is given by 

x^K, p:, A) = x: + £ d\' p:, ao) 

pli^:, p:, A) = p: - £ d\' p:, ao) . (12) 

As we will discuss in Section IV, this choice is not the most appropriate initial trajectory 
for the application to lensing. Instead, we will take a sequence of constants, (X",P^), for 
different ranges in A as the zeroth iterate. These constants will be selected from values 
predicted by the thin lens approximation. 

The nth iterate is given in the same way: X" and Pa are replaced by X^_]^(X^, P°, A') 
and P"~^(X", P°, A') on the right hand side of Eqs.(p!OD after the derivative is taken. In this 
way, the right hand sides of Eqs. (p!OD are always functions of A' and initial values, (X^, P°), 
and can be integrated immediately. The nth iterate approximate solution to Hamilton's 
equations in the original phase space variables is obtained by substituting X^ and into 
the canonical transformation, Eqs.(|^. 

In the limit that n oo, the iterative solution, (X^, P"), always converges to the exact 
solution for small A values if the initial trajectory is chosen as constant values. However, even 
in these converging cases, at any finite n, the iterative solution may differ greatly from the 
exact solution at large A values. For example, Goldstein 0] considers a harmonic oscillator 
perturbation, 

AH = -muj'^x'^, 
2 

on a free Hamiltonian in one dimension, 

2m 

Using particular initial conditions, he finds that the 2nd iterate solution in the original 
variables given by 
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X2 = {^t~ — - + — - , 

muj \ 61 5! / 

P2 = Po (^1 - ^ + ^ J , 

As n ^ oo, the limit converges to 

Po . . . 

a; — > smoji p ^ po cos out. 

muj 

However, at any finite n, the solution runs away to infinity, quickly separating from the 
exact result. This is because the finite n iterative solution is simply the first few terms in a 
Taylor series expansion of the exact result. 

To develop a method which gives reasonable results at finite tt., one may attempt to adjust 
the standard iterative method using general knowledge of the features of the exact solution. 
Such corrections are required in the application of the variation of constants method to 
gravitational lensing when accurate results at low n are desired. These corrections, which 
will be discussed in Section IV. C, involve the consecutive use of a series of unperturbed 
orbits. 

III. APPLICATION TO GRAVITATIONAL LENSING 

The iterative approach to gravitational lensing will be an approximation to the exact ap- 
proach and will remain within its overall perspective and framework In Section III. A, 
we outline the key elements of the exact approach, while the iterative approach will be 
developed in Section III.B. We show how to turn null geodesies into lens equations in Sec- 
tion III.C. For clarity, we leave any corrections of the iterative method at low n to the 
discussion in the final subsection of Section IV. 

A. Exact Gravitational Lensing 

In the exact approach, one considers a space-time, (M, containing various lenses 
(travelling on specified world lines) whose properties are hidden in the metric g""^, and an 
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observer moving on a world line described by 

x^ = X:{t). (13) 

The past light-cone of this observer is generated by null geodesies originating at X^(r), and 
may be described parametrically as 

x" = x%T,e,(j),x) 

or 

t = T(r,^,0,A) (14) 
x' = X\T,e,(f),X) i = l, 2, 3. (15) 

In these equations, the proper time, r, labels the points along the world line of the observer, 
and A represents an affine parameter along the null geodesies of the observer's past light- 
cone. The angles, {6, (p), represent the sphere of null directions seen by the observer, or the 
observer's celestial sphere. 

The time equation, Eq. ([l^) , describes the coordinate emission time of a light ray emitted 
by a source at a;* arriving at the observer at time r. This equation, often referred to as the 
"time of arrival equation," is important in astrophysics, among other reasons, because it can 
be used, in conjunction with observations and a lens model, to place bounds on the Hubble 
constant, see for example f^J^l- 

The three space-like equations, Eqs.(p!5D, can be considered as a generalized, exact form of 
the gravitational lens equations, because they map observed directions into source positions. 
By giving a relationship between the affine parameter and some observable distance scale - 
possibly the redshift distance, i? - one can express the generalized, exact lens equations 
as 

x' = A:\t,9,(I),R) i = l, 2, 3, (16) 

which is similar the traditional form of lens equations, where the (r, 6, 0, R) are all ob- 
servable quantities. In ||l|J^, Eq. ([T6|) is refered to as the exact lens equation. This step is 
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necessary for comparison wilowth astrophysical observations because the affine parameter 
has no measurable meaning, but will not be required in our discussion. 

As mentioned, the general perspective taken in exact gravitational lensing is that the 
lens properties are built into the form of the metric, g"'^. This means that there can be no 
reference to any background space-time or any quantity "in the absence of the lens." To 
find the explicit form of the lens equations and time of arrival equation, one must solve the 
null geodesic equations associated with g"''^, giving initial conditions for past directed rays 
on the observer's world line. 

Therefore, we begin with the full geodesic equations: four second order coupled differ- 
ential equations for the four base space coordinates, x". These equations can alternatively 
be viewed as a set of eight coupled first order differential equations on the cotangent bun- 
dle, {x"',Pa), which are derived from a Hamiltonian. To make them into the null geodesic 
equations, initial conditions must be chosen so that the Hamiltonian vanishes. 

The Hamiltonian for geodesic orbits is 

H = g''\x'^)paPb, (17) 
with associated Hamilton's equations: 

Pa = -g^^^aPbPc- (18) 

Solutions to Hamilton's equations have the form 

^" = 5"«,p:,a) 

Pa^Pa{xlpl,\), (19) 

where is an initial position representing the location of the observer and p° gives the initial 
direction of the geodesic. The coordinates, x", represent the space-time position along the 
geodesic, and Pa is the cotangent vector to this geodesic. 
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Null geodesies correspond to the geodesies whose value of the Hamiltonian is zero. Since 
the Hamiltonian does not depend on the parameter A, its value will be preserved along the 
geodesies, and can be set to zero by constraining the initial conditions: 

Hix:,p:) = g^\x:)p:p: = Q. (20) 

Equation ( pOD is quadratic in and can always be solved for p^ = ]5o(x",p°). Substituting 
this into Eqs.(|T9|) gives 

x'^ = x'^{x^,,p",{xlp^),plX) = x'^ixlplX), (21) 

which are null geodesies. 



B. Iterative Method for Null Geodesies 

In general, it is impossible to solve Hamilton's equations, Eqs. ([T8D , in closed form. How- 
ever, if one is able to consider the metric of physical interest, to be close to some metric, 
(7°*, for which one can solve Hamilton's equations, then the method of variation of constants 
gives a way to systematically approximate the null geodesies. 

Formally, we suppose that we have two space-times, (M, and (M, where 

g-b = gab + f^ab^ (22) 

and h"'^ is "small." For the sake of discussion, we will suppose that (M, g'^^) represents 
either Minkowski space-time or a Robertson- Walker geometry, since these are the usual 
background space-times in gravitational lensing. We form the associated Hamiltonians in 
both space-times: 

Hoix^,Pa)=gfix'')PaPb (23) 
H{x'',Pa) = g''\x'')paPb = gf {x'')PaPb + h'^' {x'')PaPb. (24) 

By finding a complete solution to the Hamilton- Jacobi equation in {M,g^^), 
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^''d^d^^^-^' ^^^^ 

we obtain a generating function of canonical transformations, Fo(x", P^, '^)- Because 
Minkowski space-time and the Robertson- Walker geometries are integrable geometries, we 
will always be able to solve Eg . (1251) and find an explicit canonical transformation in the form 
of Eqs.(0): 

= x"(X",P„,A) (26) 

Pa=Pa{X%Pa,X). (27) 



Using the explicit form of the canonical transformation, Eqs.(^, we transform the 



Hamiltonian in the perturbed space, Eq. (|2^) , obtaining 

dF 



dX 

(/^'^^(F(X^P„A))) p,(X^P„A) p,(X^P„A). (28) 



Hamilton's equations of motion for X" and Pa are 



x'^(x^ p„ A) = 2 {h'^ (^«(x^ p„ A))) p,(x^ p„ a) (^(^^ Pa, a)) 



+ ( Pa, A)) ) P„ A) p,(X^ P„ A) 



5P„ 

^S-^IX^P^A) (29) 
P,(X^ P„ A) = - (—^ (x"(X^ P„ A)) J p,{X^, Pa, X ) p,(X^ Pa, A) 

(-2/1^'= (x'^(x^ p„ A))) p,(x^ p„ A) Pa, x)^ 

= TaiX^,Pa,X). (30) 

We now use the iterative procedure of Section II to solve these equations. For the 
moment, we will choose the zeroth iterate as X" and Pa taking constant values, (X",P°), 
under the integral. The nth iterate will be given by 
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X„"(X:, P°, A) =X^+ dX' ^\XU, P:-\ A') (31) 

p:{x:, p:, a) = p:+ [ dx' t,{x:_„ p:~\ ao, (32) 

where and -P^ are functions of initial conditions and A'. 

When the nth iterate is substituted into Eg . (^61) , the result is an approximate solution 
to the geodesic equations of (M, with initial conditions (X",P°): 

= F(x„"(x:, p:, A), p:(x:, p:, a), a). (33) 

These approximate geodesies may not be approximately null. In analogy to the exact case, 
we should require that 

H{X:, P:) = 0, (34) 

which we can solve for Pq as a function of the other initial conditions, 

Po = P°(X„^i^°). (35) 

Substituting P^{X^, P") into Eq.(|33D gives approximate null geodesies: 

= x^iX^iX:, P°, A), P:(X:, P°, X),X) = x'^iX^, P°, A). (36) 

It is important to realize that, while the iteration method introduces a "background 
space-time," {M,g^^), the null geodesies one seeks are the null geodesies of the physical 
space-time, {M,g"-^). The spirit of exact gravitational lensing is that no reference to the 
background should be made. In the iterative approach, a background is introduced only to 
facilitate finding the approximate geodesies of (M, g"-^) and plays no further role. All physical 
quantities are computed in the physical space-time using g"-^, not (7"^; the null condition is 
fixed by solving Eq.(|3^) for Pq, in which the physical metric appears. Likewise, angles 
are to be computed from the inner product of two vectors which lie in the tangent space 
of (M, g"''^) . Hence, the iterative approach remains within the spirit of the exact approach 
because it does not refer in any physical way to a background space-time. 
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The nth iterate, approximate null geodesies can thus be used as the basis for an ap- 
proximate gravitational lensing theory under the spirit of the exact approach. For these 
purposes, the X" in Eq.(p6|) must be chosen as the space-time position representing the 
location of an observer, and the P° are used to parametrize the directions which an observer 
can see. Sources are then located at x", along the past directed null geodesies described by 
Eg . (p6l) . In the next subsection, we describe the procedure for turning (either an exact or 
an approximate) solution to the null geodesic equations into the form of a time of arrival 
equation and lens equations given by Eqs.(|lJ, |T5|). 



C. Exact and Approximate Lens Equations 

In this subsection, we describe how either an exact or an approximate solution to the 
null geodesic equations in an arbitrary space-time gives rise to gravitational lens equations. 
The null geodesies, exactly or approximately, are described by 

x^ = x^{xlplX) (37) 

Pa=Pa(x^,p-,A), (38) 

where the is an initial space-time location, the three p° describe the spatial direction 
of the null geodesic at that point, and A is an affine parameter along the null geodesic, 
where A = corresponds to the initial position. The null tangent vectors to these geodesies, 
pa _ go-bp^ _ j^a^ taken as past directed. The scaling of A has not yet been chosen. 



To turn Eqs.(|37|) into a time of arrival equation and lens equations in the form of Eqs.(14, 



T5|) , we introduce an observer by taking a time-like curve parametrized by the proper time. 



x^ = X:iT), (39) 

as a one parameter family of initial positions in Eqs.(|T7D, or = X"(r), for the past directed 
null geodesies. 
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The tangent vector to this world hne, 

dr 

can be used to scale the p° (and consequently A). Because the magnitude of the p° is 
irrelevant, we can fix its length by requiring that 

v^pl = -1. (40) 

This implies that, at any point along the observer's world-line, the three p" are functions of 
only two parameters. These two parameters can be taken as two independent "observation 
angles," (6',0). Thus, the null cotangent vector at X^{t) can be expressed as 

p:=p:ieA). (41) 

By replacing with X"(r) and p" with Eq.(^), in Eq. ([37|) , one obtains a time of arrival 
equation and lens equations in the form of Eqs.(|l^, p!5|): 

t = x\x:{T),pi{e,<j>),x) 
x' = x\x^{T),pi{e,<f>),x) 

= X\T,e,(j),X) i = l, 2, 3. (42) 

As an example, we consider an observer moving along a world line given by X^{t) = 
(r, 0, 0, +Zo) in a linearized Schwarzschild space-time. From Eg . (^OD , we see that Pq{x'^, pf) = 
— 1 for a past directed null geodesic. This normalization can be obtained in the previous 
section by dividing all components of Pa by the constant value of Pq given in Eq. ([35|) . 

Since p° is null, with Pq = — 1, its spatial part will have the length: 

\Pi\ = ^-g'^Pip"- = + (43) 

Thus, we can write the spatial part of the null covector at the observer as 
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where Ui{9, 4>) is a unit vector parametrized by two observation angles, {9, 4>). 

Prom the spherical symmetry of Schwarzschild space-time, the lens, source, and observer 
can be taken as lying in the same spatial plane, and any ray with an initial direction in 
this plane will remain in the plane. Fixing this plane as the x-z plane fixes one of the two 
observation angles, say = 0. An observer on the +z axis will observe the lens located at 
the origin in the direction of the unit spatial vector. 



-9ijKxi yi-© 

In terms of Ui and n*, the remaining observation angle, 9, is defined by 

cos9 = Uin\ (44) 
Since is proportional to 5*^, we have 



/I 2m 
Uz = —\ 1 1 — — r cost^. 

Zn 



Using —g^^ Ui Uj = \u\ = 1, Uy = 0, and the linearized metric given in Section IV.B, one 
finds the final component on the momentum. 



Ux — \l I — { I — T— T I cos^ 9. 



We note that the observation angle is determined relative to the physical metric of the 
space-time, . This is consistent with the notion that the iterative and exact methods 
should not make reference to a background space-time. 



IV. THE SCHWARZSCHILD LENS 



To demonstrate the iterative method, we apply this procedure to the case of the 
Schwarzschild lens, considering the Schwarzschild metric as a perturbation of Minkowski 
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space-time. Throughout this section we will work in (t, x, y, z) coordinates which, while 
awkward for the Schwarzschild metric, facilitate the iterative procedure and comparisons 
with the thin lens approach. In the iterative calculation, we work with the linearized version 
of the Schwarzschild metric, since there are no significant differences for the large impact 
parameters usually discussed. 

The section will be divided into three subsections. First, we outline the thin lens ap- 
proximation and its application to the Schwarzschild lens. The second subsection develops 
the application of the iterative method to the Schwarzschild lens. Finally, we discuss the 
choice of the zeroth iterate when one wishes to use low n values. This subsection makes use 
of results from both of the first two subsections. 

A. The Thin Lens Approximation 

Here, we briefly outline the standard thin lens approximation. Most theoretical de- 
tails about the thin lens approximation can be found in the book by Falco, Schneider and 
Ehlers §. 

The key assumption in the thin lens approximation is that the null geodesies in the 
space-time of interest travel along the null geodesic path of the background space except 
at isolated points near the lens where sharp bending occurs. For a Schwarzschild lens in a 
Minkowski background, the geodesic travels in rectilinear motion from the observer, O, to a 
point in the lens plane, I, and then from I to the source location at O, as is shown in Fig. |l|. 

The effect of the lens is twofold. First, the trajectory is instantaneously bent at I by 
a bending angle computed as the angle between the two asymptotic paths of the true null 
geodesic. For a Schwarzschild lens, the bending angle is simply a = AG m / {(? To), where To 
is the radius of closest approach, m is the mass, and the physical factors G, and c have been 
restored. Second, the lens affects the time of arrival by introducing a gravitational time 
delay relative to the time required to traverse the same path in the Minkowski background. 

The lens equation (of the thin lens approximation) is the map from the space of observed 



16 



light-ray directions to the source positions (on the source plane). Often the map is not 
invertible, as there could be multiple images for a given point source. In practice, one 
develops a model for the lens, usually in terms of a potential for the mass distribution, and 
from this model determines the amount of bending a null geodesic will encounter along its 
trajectory from past null infinity to future null infinity. This bending angle, a, in Fig. is 
amount of bending of the ray at the lens plane. An angle (3 locates the true source location 
in the background space-time. 

For a single Schwarzschild lens, one determines from Euclidean geometry that the lens 
equation is 

AGmDis , . 

^^'-^^^DMi' '''' 

where 9 is the image observation angle, and Di, Ds, and Dis represent the distances between 
the lens and observer, the source and observer, and the lens and source, respectively. Gener- 
ically, there will be two values of 6 for a given value of (3, so that there will be two distinct 
paths joining the source and observer, with one passing on each side of the lens. 

In the thin lens approximation, the time of arrival is determined by the coordinate time 
which elapses along the path taken in the thin lens approximation. Although the path, 
up to a sharp bending in the lens plane, is determined by the background space-time, the 
time of arrival is computed relative to the (usually linearized) physical metric including the 
lens. The time of arrival for a ray emitted from the source at time ts and passing close to a 
Schwarzschild lens is 

1 r f 2Gm\ „ , , 

where r{l) represents the Euclidean distance from the lens to a point along the trajectory. 
This distance is expressed as a function of the Euclidean distance along the trajectory, I, 
and the integral is taken along the rectilinear path from O to I to S. A time delay is the 
difference between the two times computed by Eq.(^) along the two paths determined by 
Eq.dl). 
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B. Iterated Schwarzschild Null Geodesies 



This subsection develops the iterative approach in the case of the single Schwarzschild 
lens along the lines of Section III.B. In this subsection, the zeroth iterate will be chosen as 
the solution to Hamilton's equations in the unperturbed metric, or {X^, P°) = constant. A 
better choice of the zeroth iterate is proposed in the next subsection. We begin with the 
linearized Schwarzschild metric, which is given by 

/ 
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(47) 



in {t, X, y, z) coordinates. Here, the signature of the metric is (+ 



5 ! 7 )■> 



and r 



\J x"^ + y"^ + z"^ . The Hamiltonians of Minkowski space-time and the linearized Schwarzschild 
space-time are 



1=1 



(48) 
(49) 



and h"'^ is 



r 

2mx^x^ 



(50) 



The complete solution to Hamilton- Jacobi equation of Ho{x"',Pa) is 



Fo{x^,PaA) = x''Pa-v"'PaPbX 



(51) 



and the canonical coordinate transformation associated with Eq.(^l|) is given by 



x'' = X'' + 2r/"^PbA = F(X^ Pa, A) 

Pa = Pa- 



(52) 
(53) 



Applying this canonical transformation to the Hamiltonian H{x"',pa) yields 

i^'(x^ p„, A) = h""" p„, A)) n, (54) 

and Hamilton's equations of motion for X°- and Pa are 

X-{X% Pa, A) = 2h-' Pa, A)) n + 9h'HxY^,Pa,^)) p^ 



dP, 



a 



= E\X\Pa,\) 

Pa{X\ Pa, A) = JhHnX^.PaA)) p^ ^ ^ ^^^a ^ p^^ (55) 

We solve Hamilton's equations by the iterative procedure, as in Section HI. The zeroth 
iterate will be constant values, X° = X^ = and Pa = Pa = Pa- The nth iterate is given 
by 

x:ix:, p:, a) = x: + f d\' e%x:^„ p:-\ ao 

p:{xi p:. A) = p:+ j dx' t,(x^„ pr\ ao. (56) 

Jo 

The approximate geodesies given by Eq. (|56D are not necessarily null; we must require that 



H(X^,P°) = 0. Using the metric, Eq.(47), and the Hamiltonian, Eq.(49), yields 



where = J 5ij XlXi. When Eqs.(^) and Eq.([57|) are substituted into the canonical 



transformation, Eqs.(52), the nth iterate, initially null geodesies are obtained. 

Finally, we note that when the initial conditions satisfy the null condition, Eq.(^), 
the resulting trajectories are approximately null in {M,g°'^). Geodesies in the unperturbed 
background with these same initial conditions will be space-like. Hence, a cone which is null 
in (M, g'^^) will be fully inside the light cone of (M, g"'^) due to the converging effect of the 
positive mass in the Schwarzschild lens. 
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C. Choice of the zeroth iterate for low n 

As mentioned at the end of Section II, the iterative method can fail to produce accurate 
results at low n and large A values regardless of whether the method converges in general. 
There are two ways in which such a failure can be anticipated in the application to lensing. 

First, when the zeroth order trajectory is chosen as one set of constants, X" = X" and 
Pa = Pa, the trajectory is a "straight line," corresponding to a solution to the equations 
of motion of the unperturbed Hamiltonian. In the Schwarzschild case, the first order tra- 
jectory is a "curve" which bends in the vicinity of the lens. These two trajectories become 
increasingly separated as A grows. For n = 1, the integrands in Eqs. (|56|) are evaluated along 
the zeroth order trajectory, and the "force" acting on the "photon" becomes increasingly 
incorrect as A grows. Hence, when the zeroth order trajectory is far from the first iterate 
trajectory, one may anticipate poor results. 

It turns out that for a single Schwarzschild lens, this particular source of error is reason- 
ably insignificant because the value of the "force" is decreasing faster than 1/ A as A grows. 
However, in a situation where there are multiple lenses lying in more than one lens plane, 
the error becomes very significant, and the force vector may point in the wrong direction 
altogether. 

A second way to anticipate troubles is to note that while the true Hamiltonian is con- 
served under a flow along the actual geodesic, it will not be conserved under an approximate 
flow. If one finds the exact solution for (X", P"), and sets the initial conditions as in Eq.(^), 
the value of the perturbed Hamiltonian, Eq . (p4|) , will remain zero along the trajectory. How- 
ever, at any finite n, the value of 

H{X:, P°, A) = H{ i^^ A), p:{X:, P^, A) ) (58) 

will drift away from zero because the nth iterate solution is not truly a null solution to the 
geodesic equations. Hence, at large A, the nth iterate solution may not closely approximate 
a null geodesic. 
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The basic reason for these problems is that, so far, we have chosen zeroth iterate in a 
non-physical way, as a straight line path corresponding a single solution to the equations of 
motion for geodesies in the unphysical space-time, (M, This solution separates from 
the first iterate shortly after passing the lens plane. 

The simplest choice of the zeroth iterate is the solution of the equations of motion of 
the unperturbed Hamiltonian, but on the other hand, the zeroth iterate should remain close 
to the first iterate trajectory. As we have discussed, the thin lens path is given by two 
connected solutions of the equations of motion in (M, g'^^) and remains close to the exact 
solution. 

Therefore, we may avoid these problems by choosing the thin lens path as the zeroth 
iterate. This is accomplished by choosing different sets of constant values, (X",P°), in 
different A ranges as described below. Making this choice insures the first iterate trajectory 
will always remain close to the zeroth iterate regardless of the number of lenses. 

From a kinematical standpoint, the entire thin lens trajectory is determined by the choice 
of X* and P°. These values determine the straight line path from the observer at O to I 
in Fig. They also determine a new set of six parameters, (X^^, pI^) which represent the 
location of the intersection of the original straight line path and the lens plane, {Ip), and 
the new momenta determining the straight line path joining I and S. 

As Schwarzschild space-time is static, the X° = T and Pq equations separate from the 
spatial equations. Due to spherical symmetry, we can assume that the spatial path of the 
geodesic lies in the x-z plane. In this case, the thin lens trajectory is determined by 



x = Xl-2 P°A 



< A < A, 



(59) 




2^3 A 



< A < A, 



(60) 



and 



X = Xl - 2 P^^A A; < A < A, 
z = Xf - 2 P^^A Az < A < A, 



source- 



source 



(61) 



(62) 
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If the initial conditions represent an observer on the z axis at X\ = 0, and = Zo, the 
lens plane will be 2; = 0, and we have 

A; 



2 Pi 



— P° 7 

^Ip - 

4 = 0- (63) 
The scaling for A and the spatial part of the Pa is determined by 

v'^Pa = -1, (64) 

as in Section III.C. Then, by choosing the new P/^ such at the angle between P/^ and the 
original P° at Xl^ is the bending angle, a = AG m/{c^ Vo), 

ij pip po 
^oi j 



COS a = — I ' , (65) 



both the i^'^ are determined. 

In summary, we may always pick the initial conditions as {Xl = 0,X^ = z^yP" = 
p°, P3 = This fixes the value of A; reaching the lens plane and Xl^. Eg . (|6^ and Eg . (|65D 
determine the P-^ in terms of {zo,pl,p°). Hence, by freely choosing initial values at the 
observer, the entire zeroth iterate path is determined. 

The first iterate solution will then be represented by the integrals, Egs.(^), where for 
< A < A; the constant values are fixed at the initial point, and for A; < A < Xsource, the 
values are reset. The initial Pq is fixed as in Eg. ([57|) at the initial point along the first leg, 
and then reset at X/^ with the new values. This insures that the value of the Hamiltonian, 
H = g""^ PaPb will remain close to zero along the trajectory. 

As an example, in Fig. ^ we plot a first iterate trajectory where the initial conditions 
are taken from the thin lens path as mentioned in this section. This figure also shows a 
thin lens trajectory with similar initial conditions in dashed lines. The mass is taken to be 
1.0 in geometrical units and the range in z is 100 units. With these parameter choices, the 
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differences between the first iterate and thin lens are greatly enhanced compared to what 
one would see for a more physical choice of distance scale. 

We see that the first iterate trajectory is a smooth path, while the thin lens path has a 
kink at the lens plane. This demonstrates one fundamental difference between the thin lens 
and iterative approaches, namely, that the iterate approach represents a continuous bending 
process while the thin lens approach reduces the action of the gravitational force to a single 
lens plane. We also note that there can be a significant difference in the distance of closest 
approach between the two methods, which leads to a difference in the bending angle and 
possible location of a source. 

The nth iterate time elapsed along the trajectory is formally given by 

tn = to + [ dX2Po(l + ^), (66) 



rn-i_ 

from Eg. (|56|) , the form of the perturbation, /i"'', and the canonical transformation. Here, 
r„_i is defined as the radial coordinate distance from the origin to the point along path of the 
n — 1 iterate. It is not difficult to compute the second iterate for the time by first computing 
the first iterate trajectories above, and using these in the right hand side of Eq.(^). In 
this case, the integration over A will also contain two segments, corresponding to the two 
segments above. 



V. DISCUSSION 

Gravitational lensing is a very active field of observational research and is becoming 
an increasingly important tool for probing the structure of lenses and the cosmological 
parameters. Each day, new observations and analysis expand the list of lensing candidates, 
providing a rich field of study. 

While the thin lens approximation has proven remarkably accurate when applied to the 
current observational data, it relies heavily on the weak field and small angle approximations. 
With the ever expanding hsts of candidates and improving observational techniques, it is 
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not difficult to imagine that the hmits of the thin lens approximation may be approached 
in the future. Virbhadra and Ellis attempt to improve upon the thin lens approximation by 
removing the small angle approximations 0, and other methods may try to "fix up" the 
thin lens approach from within its overall framework. An another approach was developed 
by Ted Pyne and Mark Birkinshaw |TyjTT|, where the null geodesic equations are solved 
pertubatively. One may also numerically solve the geodesic equations, as in the "direct 
integration method" of Tomita, et al ||l2| , p!3| . 

The iterative method provides a different approach. Although it works within the frame- 
work of exact gravitational lensing, the iterative approach is reasonably easy to apply to 
a wide range of lensing scenarios. In addition, when the thin lens trajectory is taken as 
the zeroth order path, the iterative method represents possible improvement upon standard 
lensing equations, with additional accuracy obtained by increasing the level of iteration. 

We are preparing a close study of the iterative and thin lens approaches in a wide variety 
of space-time settings, including multiple lenses in multiple lens planes. Of particular interest 
are the time delays in the multiple lens plane configurations. It is entirely conceivable 
that the thin lens approximation may be accurate beyond any observational limits, and no 
noticeable differences may be found with the iterative approach. However, we take courage 
in the fact that the history of gravitational lensing is ffiled with false prophecies of the 
unlikely nature of detection. 
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FIGURES 



o 




FIG. 1. The general picture used in the thin lens approximation. A past directed geodesic 
emitted at the observer O travels rectilinearly to I, is bent by the bending angle, a, and then 
travels to a source at S. Distances Di, Dg, and Dig separate the lens and observer, source and 
observer, and lens and source, while /? locates the source and 9 is an observation angle. 
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FIG. 2. The spatial path of the thin lens and first iterate approximate null geodesic when the 
first iterate takes the thin lens approximation as its zeroth iterate. 
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